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Abstract 

We prove local existence and uniqueness of static spherically symmetric solutions of the Einstein- 
Yang-Mills equations for any action of the rotation group (or SU(2)) by automorphisms of a principal 
bundle over space-time whose structure group is a compact semisimple Lie group G. These actions 
are characterized by a vector in the Cartan subalgebra of g and are called regular if the vector lies 
in the interior of a Weyl chamber. In the irregular cases (the majority for larger gauge groups) the 
boundary value problem that results for possible asymptotically flat soliton or black hole solutions is 
more complicated than in the previously discussed regular cases. In particular, there is no longer a gauge 
choice possible in general so that the Yang-Mills potential can be given by just real-valued functions. 
We prove the local existence of regular solutions near the singularities of the system at the center, the 
black hole horizon, and at infinity, establish the parameters that characterize these local solutions, and 
discuss the set of possible actions and the numerical methods necessary to search for global solutions. 
That some special global solutions exist is easily derived from the fact that su(2) is a subalgebra of any 
compact semisimple Lie algebra. But the set of less trivial global solutions remains to be explored. 
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1 Introduction 

The classical interaction between gravitational and Yang-Mills fields is described by a complicated highly 
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nonlinear field equations which, even when reduced to a system of ordinary differential equations in the static 
spherically symmetric case, leads to many interesting mathematical and physical problems. Physically these 
solutions have shown that equilibrium configurations of black holes can be much more complicated than had 
previously been thought since mass, charge and angular momentum are clearly not enough to characterize 
them. There is even numerical evidence for the existence of non spherical axisymmetric static black holes JHJ . 
Mathematically, an analysis of the solution space of the static spherically symmetric equations requires an 
interesting combination of geometrical, algebraic, analytic and numerical techniques. The global solutions 
of SU (2)-EYM equations have been extensively studied analytically |l6|, p~9|-[2T]| . For a fairly comprehensive 
summary of the substantial literature on the subject we refer to the review article |22]. Almost all these 
investigations, however, have only studied the gauge groups SU(2) and occasionally SU(n) for n > 2, and 
only for the most obvious ansatz for a spherically symmetric gauge field. 

But spherical symmetry for Yang-Mills fields is more complicated to define than for tensor fields on a manifold 
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because there is no unique way to lift an isometry on space-time to the bundle space. The only natural way 
to define spherical symmetry of a Yang-Mills field is to require that it be invariant under an action of the 
rotation group by automorphisms on the principal bundle whose structure group is the gauge group G. A 
conjugacy class of such automorphisms is characterized by a generator Ao which is an element of a Cartan 
subalgebra [) of the complexified Lie algebra g of G || . If one restricts consideration to fields which are 
bounded at the center or, in the presence of black hole fields, to those for which the Yang-Mills-curvature 
falls off sufficiently fast at infinity then A must be a defining vector of an s[(2)-subalgebra of g. 

One of these classes of actions of the symmetry group is somewhat distinguished. It corresponds to a principal 
defining vector in Dynkin's terminology and we will call it a principal action. Almost all work for larger 
gauge groups has been done for this case p^ , |l3| ,pq,p!7|]. A slightly bigger class of actions, which we call 
regular, consists of those for which the defining vector lies in the interior of a Weyl chamber. For those, 
for example, Brodbeck and Straumann ^,0] proved that all bounded static asymptotically flat solutions are 
unstable against time dependent perturbations. 

Of course, one is interested in global solutions of the boundary value problem that results from demanding 
boundedness at the singularities of the differential equations at the center or the horizon and at infinity. 
This global existence has long been established for G = SU(2) by Smoller et al. [^0[|l[ (see also ||). It is 
easy to "imbed" these solutions into the set of solutions of the theory with arbitrary compact semi-simple 
G since the latter always has SU (2) subgroups. The problem is thus not to prove that global solutions exist 
but to explore the global solution space and, hopefully, characterize different types of solutions, for example, 
by their behavior near the center or near infinity. 

In p8| we have therefore considered and solved the local existence problem for bounded solutions near the 
singularities for the regular symmetry group actions. We have identified the set of "initial conditions" that 
must be given at r = or r = oo to guarantee the local existence and uniqueness of a bounded solution. For 
arbitrary gauge groups this required an fairly intricate application of the sl(2, C) representation theory on 
the (complexified) Lie algebra g of G. 

The purpose of the present paper is simply to extend these results to the case where the symmetry group 
action need not be regular. It turns out that the situation is qualitatively not very different. There are 
still similar algebraic equations that restrict the possible initial data. But the number of functions that will 
characterize the gauge potential is no longer just the rank of the Lie algebra g or of one of its subalgebras. 
Even the set of reduced field equations can no longer be determined simply from the structure of the Cartan 
subalgebra (one needs access to all the Lie brackets of g). Moreover, a simple gauge choice that allows one 
to describe the potential in terms of rank(g) real functions is no longer available. Complex functions must 
be allowed which increases substantially the number of parameters to be determined in a numerical solution 
of the boundary value problem. 

In section |2| we review definitions and previous results mostly from |l8| and then describe in section ||| how 
the field equations can be handled computationally. In section || we give some examples of the possible 
symmetry group actions for low dimensional gauge groups. The methods of section ^ do not lend themselves 
easily to a general existence proof. This needs to be done differently. In section || we state and prove some 
algebraic lemmas that needed in section |6| for the local existence and uniqueness proofs. We conclude by 
giving a preliminary example of a numerical irregular solution showing that imaginary parts of the functions 
w a (r) develop even if the values w a (0) or w a (oo) are all real. We have not yet found a (nontrivial) global 
asymptotically flat numerical solution. 

2 Yang-Mills potentials and field equations 

Let P be a principal bundle with a compact semi-simple structure group G over a static spherically sym- 
metric space-time manifold. For simplicity we consider only actions of the group SU (2) by principal bundle 
automorphisms on P that project onto the action of 5*0(3) on space-time which defines the spherical sym- 
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Equivalence classes of these spherically symmetric G-bundles are in one-to-one correspondence to conjugacy 
classes of homomorphisms of the isotropy subgroup, U(l) in this case, into G. The latter, in turn, are 
given by their generator A3, the image of the basis vector T3 of su(2) (where {Ti} is a standard basis with 
[t»»7j] = e ijTk), lying in an integral lattice I of a Cartan subalgebra f) of the Lie algebra go of G. This 
vector A3, when nontrivial, then characterizes up to conjugacy an su(2) subalgebra. (See, for example, Q; 
we follow the notation of this text and also of fllPf.) 

It is convenient to pass to the complexified Lie algebra g of go and define 

A := 2*A 3 . (2.1) 

We now regard go as a compact real form of g which defines the conjugation c on g (a Lie algebra automor- 
phism satisfying c o c = I ). Then we can write g = go + igo, c(X + iY) = X — iY if X, Y 6 go, and also 
fj = fjo + ifyo where f)o is a Cartan subalgebra of go- Moreover, 

A e it) , c(A ) = -A . (2.2) 

We will use the following notation from Lie algebra theory (following [Q, ||, Q): Ad 9 : go — * go V g € G 
is the adjoint action of the Lie group G on its (real) Lie algebra go while ad : g — > g[(g) is defined by 
ad(X)(Y) := [X, Y]. Define also the centralizer of X in g by 



{Feg| [X,Y]=0] (2.3) 



and write g x for the corresponding centralizer of the real Lie algebra. 

Wang's theorem (jl4],^3| on connections that are invariant under actions transitive on the base manifold has 
been adapted to spherically symmetric space-time manifolds by Brodbeck and Straumann ||. They show 
that in a Schwarzschild type coordinate system (t, r, 6, 4>) and the metric 

g = -NS 2 dt 2 + N^dr 2 + r 2 {d0 2 + sin 2 Odcf 2 ) (2.4) 

a gauge can always be chosen such that the go-valued Yang-Mills-connection form is locally given by A — 
A + A where 

A = N(t, r)S(t, r)A(t, r)dt + B(t, r)dr (2.5) 

is an Ad(A3)-invariant 1-form (i.e. with values in Qq 3 ) on the quotient space parametrized by the r and t 
coordinates and 

A = AxdO + (A 2 sine + A 3 cos 0)# . (2.6) 
Here A3 is the constant isotropy generator as above and Ai and A2 are functions of r and t that satisfy 

[A 2s A 3 ]=Ai and [A3,Ai] = A 2 . (2.7) 

In this paper we will only consider the static magnetic case for which Ai and A 2 as well as N and S are 

functions of r only and the 'electric' part A of the potential vanishes. 

With 

A±:=T-Ai-iA 2 (2.8) 



1 This ignores some interesting effects due to the fact that 50(3) is not simply connected. For an analysis of SO(3) actions 
on St/ (2) bundles see 0. 
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equations (2/7) become 

[A ,A ± ]=±2A± (2.9) 

and the A±(r) satisfy the reality condition 

A T = -c(A±) (2.10) 

With this choice for the gauge potential, the Einstein- Yang-Mills (EYM) equations can be written as j| 

m' = NG + r~ 2 P, (2.11) 

S- 1 S , = 2r- 1 G, (2.12) 

r 2 NA+ + 2(m — r~ 1 P)A' + + J 7 = 0, (2.13) 

[A' + ,A_] + [A'_,A + ]=0 (2.14) 



where ' := d/dr and 



N=:l-— , G:=±{A' + ,A'_), P:=-\{F,F), (2.15) 



F:=i(A -[A +) A_]), (2.16) 
F :=-i[F,A+]. (2.17) 

Here ( , ) is an invariant inner product on g. It is determined up to a factor on each simple component of a 
semi-simple g and induces a norm |.| on (the Euclidean) [) and therefore its dual. We choose these factors 
so that ( , ) restricts to a negative definite inner product on go. 

For several purposes, in particular for numerical solutions, equations ( |2.11 )-( [2T3 ) are best replaced by an 
equivalent system that regularizes the almost singularity when N is close to zero. These equations, introduced 
by Breitenlohner, Forgacs and Maison g in the SU(2) case, take the form 

r = rAf, AT = Af(K - AT) - 29, it = 1 - K? + 29, S = (K - Af)S (2.18) 

A' + = rU+ and U+ = -{K - N)U+ - Tjr (2.19) 

where 

Af:=VN, S:=MS, K := ^ (l + N + 2G - 2P/r 2 ) , U+:=NA' + , Q := \\\U + \\ 2 (2.20) 

and the dot denotes a r-derivative. Note that the r variable behaves somewhat like the logarithm of r near 
and for r — > oo (where Af — * 1). 

For later use, we introduce a non-degenerate Hermitian inner product g by (X\Y) := — (c(X),Y) for all 
X,Y in g. Then ( | ) restricts to a real positive definite inner product on go- From the invariance properties 
of ( , ) it follows that ( | ) satisfies 



(X\Y) = (Y\X), (c(X)\c(Y)) = {X\Y), and ( [X, c{Y)]\Z) = {X\[Y, Z] ) 

for all X, Y, Z G g. Treating g as a M-linear space by restricting scalar multiplication to multiplication by 
reals, we can introduce a positive definite inner product (( )) : g x g ^ 1 on g defined by 

({X\Y}} :=Re{X\Y) V X,Y e g . (2.21) 
Let || |j denote the norm induced on g by (( | )), i.e. 



|LY|| := y/((X\X)) VIeg. (2.22) 
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From the above properties satisfied by ( | ), it straighforward to verify that (( | )) satisfies 

((X\Y)) = ((Y\X)), ((c(X)\c(Y))) = {{X\Y)), and (( [X, c{Y)]\Z » = {{ X\[Y, Z] » (2.23) 
for all X,Y,Z £ q. 

Using the norm (2.22), equations (2.15) can be written as 



G = i||A' + || 2 



and F=i||F|| 2 



(2.24) 



which shows that G > and P > 0. Energy density, radial and tangential pressure are then given by 
Aire = r~ 2 (NG + r- 2 P), Airp r = r~ 2 (NG - r" 2 P), inpe = r~ 4 P. 



(2.25) 



The main result of this paper is that the EYM equations ( 2.11 )-( 2.14 ) admit local bounded solutions in the 
neighborhood of the origin r = 0, a black hole horizon r = rjj > 0, and as r — > oo. To prove this local 
existence to the EYM equations (2.11)-(2.14), we proceed in three steps. 

1. First we prove the existence of local solutions {A + (r),ni(r)} to the equations (2.11) and ( 2.13| ). 

2. Then we determine which solutions from step 1 satisfy equation (2.14). 



3. Finally, equation (2.12) can be integrated for all solutions from step 2 to obtain the metric function 
S(r). 

Our method for carrying out the first step will be to prove that there exists a change of variables so that 
the field equations (2.11) and ( [2.13] ) can be put into a form to which the following (slight generalization of 
a) theorem by Breitcnlohncr, Forgacs and Maison || applies. 

Theorem 1. The system of differential equations 

diij 

t^ = t^Mt,u,v) i = l, 



, m 



t-r- 

dt 



-hj(u)vj + t Vi gj(t, u, v) 



. n 



(2.26) 
(2.27) 



where fj,i, Vj are integers greater than 1, and gj analytic functions in a neighborhood of (0, Co,0) G 
M 1+m+Tl , and hj : R m — > K functions, positive in a neighborhood of Cg € R m , has a unique analytic solution 
1 1 — ^ (ui(t),Vj(t)) such that 



Ui(t) =Ci + 0{V H ) and Vj (t) = 0{t v ') 



(2.28) 



for \t\ < R for some R > if |c — Co| is small enough. Moreover, the solution depends analytically on the 
parameters Cj. 

The next lemma shows that if {A + (r),m(r)} is a solution to the field equations (2.11) and ( 2.13j ) then the 
quantity [A' + , A_] + [A'_, A+] satisfies a first order linear differential equation. This unexpected result is what 
allows us to carry out step 2 and thereby construct local solutions. 

Lemma 1. If{A + (r),m(r)} is a solution to the field equations ( 2.1l| ) and ( |2.13| ) then^(r) := [A+ (r) , A'_ (r)] + 
[A_(r), A' + (r)] satisfies the differential equation 7' = — 2(r 2 N)~ 1 (m — r _1 P) 7 . 

Proof. Differentiating 7 yields 
7 / = [A + ,A^ + [A_ ) A'|] 



2 N 



P 7. 



r 2 N 



([A_,[F,A + ]] + [A + ,[A_,F]]) 



by (gTlj) and ( p,13| ) 



while [A_, [F, A + }] + [A + , [A_,F]] = by (p.9|), ( |2.16| ), and the Jacobi identity. Combining the these two 
results proves the lemma. □ 
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3 Solving the field equations computationally 



Here we describe without proofs the more practical aspects of solving the field equations. That all these 
constructions work will follow from the proofs given in sections ^ and ^. 

We consider only situations where the EYM field is nonsingular at the center and/or the gravitational and 
Yang-Mills field fall off rapidly at infinity. From the expressions for the physical quantities ( 2.2 5| ) and (2.24) 



it then follows that F vanishes there so that by ( 2.1 6| ) {Ao, A+, A_} form a triple defining an st(2) (or A 



subalgebra of g. This restricts the (constant) Ao vector to be a defining vector of such a subalgebra which 
have been classified by Dynkin ||. Alternatively, in the terminology of ||, the set {Ao, Q + , Q-}, where we 
define Q± to be the limiting values of A± at r = or oo, is called a standard triple defining a nilpotent orbit, 
A is the neutral and Q + the nilpositive element. 

It is known |J that there is always an automorphism of g that maps A onto the fundamental Weyl chamber 
so that its characteristic (or weighted Dynkin diagram) x = (Xii--- iXt) '■= (cti(Ao),... , af(Ao)) with 
respect to a chosen Cartan subalgebra and its dual basis {a±, . . . , ag}, satisfies \k G {0, 1, 2} V k = 1, . . . ,£. 
We call the symmetry group action and the vector Ao regular if Ao lies in the interior of the Weyl chamber 
or, equivalently, if all Xk are positive. The action and Ao are called principal if Xk — 2 Vfc. 

Since for semi-simple gauge groups all constructions are easily decomposed into those of the simple factors 
a classification need only be done for the simple Lie groups. It turns out that only for the Ag (or sl(£ + 1)) 
series of the classical Lie algebras and then only for even £ there are regular actions other than the principal 
one. (See Theorem 2.) Moreover, in those cases the gauge field corresponds to one of a direct sum of 
two lower-dimensional simple Lie algebras of type Ag. It is for that reason that we could confine ourselves 
in Jl8| to the principal case when studying the local existence problem of solutions for regular symmetry 
group actions. 



In general, given a fixed defining vector Ao, the Yang-Mills field will be fully determined, in view of ( 2.10 ), 
if A + is known as a function of r. Condition (|1|) implies that A + must lie in the vector subspace V2 of g 
where we define, more generally, the eigenspaces of ad(Ao) by 

V n := {X G 3 I A .X = nX} neZ, 
Here, and in the following, the adjoint action of spanc{Ao, f2+, fl-} = st(2, C) on g, is denoted by a dot, 

X.Y := &d(X)(Y) Vie span c {A , Q+, Y G g. 

In terms of a Chevalley-Weyl basis (see jl(| or jl8| for the notation) {hj := h aj , e Q , e_ Q |j = 1,... ,£;a € 
R + }, where R + is the set of positive roots of g with respect to the root basis {ot\, . . . , ag}, we then have 

A +0) = ™«( r ) e " e v 2 and A -0) = w a {r)e^ a S V-2 (3.1) 



where 



S Ao :={aeR\a(A ) = 2}. (3.2) 



In the regular case the Stiefel set 5a is necessarily a n-system ||, i.e. forms the basis of a root system 
generating a subalgebra Q\ of g (sec B). Since for base root vectors a,/3 also [e^e^] 7^ only if f3 = —a 



and [e a ,e_ a ] = h a it follows then from (2.16) that F lies in the Cartan subalgebra of gA . Substituting 



fl3T| ) int0 Q2-14j ) leads to 

w' a wp - Waw'p = Va, (3 g Sa (3.3) 

so that all the phases of the complex w a (r) are constant. 



G 



In the general case, Sa will not be linearly independent and the set {[e^e-^] | a, S S\ } will not lie in 
the Cartan subalgebra f) of g. There is no simple gaug e transformation now that will allow us to choose the 
real, although, as follows from Lemma [j], if ( |2 .14 ) is satisfied at one regular value of r it will be in a 



whole interval 



The Yang-Mills field F no longer takes all its values in the Cartan subalgebra. It is instead given by 

- i w a (r)wp(r)[e a ,e-p] . (3.4) 




However, the term T appearing in the Yang-Mills equation (2.13) lies in V2 since [X, [Y, Z]] 6 V2 whenever 

two of the X,Y,Z lie in V2 and one in V—2- 

For computational purposes ( 2.13| ) can then be written 



where 



r 2 Nw^ + 2(m-r- 1 P)w' Q + f a =0 (3.5) 



fa ■= W a + - ^ (J'a^ySWfiWyWg (3.6) 



and 



- p =^ll A o|| 2 - ^ \a\ 2 \w a \ 2 -^ ^ \a\ 2 ^ap^5WaWf}W 1 WS- (3.7) 
oGSaq a,/3,7,<5eSA 

(Here the Jacobi identity and the invariance of the inner product ( , ) were used.) Thus the general structure 
of the Lie algebra q enters the equations only via the quantity fiap-ys which is defined by 

yp n ,SeS Ao . (3.8) 

To start off the numerical integration of a bounded solution near one of the singular points we need to sum 
a power series at a point nearby. This is again similar to the method used in the regular case, but we must 
allow for possibly complex functions w a (r). This has most of all the unpleasant effect that the values that 
the w a can take at the center and at infinity form no longer just a finite set with a few signs to be chosen 
arbitrarily, but an Af-dimensional real variety in the space C M (where M denotes the number of elements 
in S\ and thus the number of functions w a needed to characterize the gauge potential A + ). Since there is 
no reason to believe that for a global solution on [0, oo) the values of A + (0) and A + (oo) should be the same 
we can, for example, choose an arbitrary fl + = A + (oo) which amounts to a global gauge choice. But then 
the value of A + at the center could be any point of this M-dimensional real variety so that the coordinates 
describing the latter need to be given as parameters as well as some of the first few power series coefficients. 



Wishing to find a local analytic solution to equations (2.11) and (|3.5|) we expand all quantities in a power 



series in r near r — 0. The lowest order terms w a fi — w a (0) are constrained by F = so that 

A = [f2 + ,ft_] and = -c(0 + ) or (3.9) 



A =:^Aihi= ^ w a ,oWi3,o[e a ,e-p] (3.10) 

i=l a,/3eS An 
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which gives, if we write h Q = 53i=i h a ,i^i and R denotes the set of all roots of the Cartan subalgebra () of 9, 
^2 h a ,i\ w a,o\ 2 = \ (i = l,... ,t) and ^ w afi wpfl[e a ,e-p]=Q. (3.11) 

aeSA q,/3gSa 

a-0eR 



We have not yet been able to determine whether equations ( 3.11 ) have solutions for all simple Lie algebras 
and all defining vectors Aq. For all low-dimensional examples we have so far considered solutions exist and 
form an M-dimensional real variety in C M . Moreover, it appears that there always exist vectors f2+ with 
only real components with respect to the basis {e Q }. 

Suppose now that f2+ has been chosen. Then equation ( [2.11 ) yields the recurrence relation 



m k = i ^P fe+ i + G k -! - 2jj n k-iGiJ (3-12) 

and (JDJ) gives 

A(A+. fc ) - k(k - l)A+, fc = b fe (fc = 2,3,...) (3.13) 
where the are complicated expressions of lower order terms and A is defined by 

A= ±ad(fi + )° (ad(f2_) + ad(f2 + )oc). (3.14) 



This A is an E-linear operator which turns out to be symmetric with respect to the inner product (2.21) 
and restricts to V2 ■ As will be shown in section ^ half of its eigenvalues are zero and the remaining ones are 
positive integers of the form k(k + 1) for fc = 1, 2, We will later need the notation 

£ := {s G N I s(s + 1) is an eigenvalue of A} (3.15) 

Moreover, to every eigenvector with nonzero eigenvalue there exists one (its multiple by i) that has eigenvalue 
zero. 

Explicitly, if we write f2+ = Y^J=ii u a,o + w tt ,o) e Qi where M := \S^ \, then A is given by the matrix 

(A)=^ (3.16) 



where 



(7 



a/3 



^ (^q5(/3 7 ) u 7,0 W (S,0 + Ma<5[/37] u 7,0^o) (3-17) 
7,<5 

b a p — ^(Ma5[/3 7 ] - ^aj(pS)) u f,O v S,() (3.18) 
7,<5 

Cap = ^(Mq7[,35] - Mq5(^7)) U 7,0 W (5,0 (3.19) 
7,<5 

da(3 = — ^ (Vaylt38]Uj,0US,0 + ^a 1 ((3S)V 1 fiVgfi) (3.20) 
7,<5 

Now let E be the matrix whose columns are eigenvectors of A such that 

AE = EA (3.21) 
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where A = diag(0 . . . OAi . . . Am) with Aj < Aj if i < j. Then we know (from section ||) that if E is of 
the form E = then also will be a matrix of eigencolumnvectors with now the correct 



correspondence between the real and imaginary parts. 

To generate the powerseries for u and v we now find from the recurrence relation ( [3. IS ) for the coefficients 
of r k that 

X fc ) = * TXk (3 - 22) 

where Xk. a — ((E T ) _1 bfe) Q /(A a — k(k — 1)) if A Q 7^ k(k — 1) or a new free parameter otherwise. 

This shows that the general solution near r = is of a similar form as in the regular case ( |18|, Theorem 
4), namely 

M 

w a (r) = w afi + Y / ^ a 0r k ^ +1 w f3 {r), a=l,...,M (3.23) 

0=1 

where k a (k a + 1) is the a-th nonzero eigenvalue of A and w a (r) are analytic functions near r = 0. The 
general solution {m(r), w a (r)\a = 1, . . . , M} is determined by M real parameters for the initial values w a (0) 
as well as another M real parameters that can be arbitrarily chosen in the coefficients ^(O) of r k,3+1 . 
The general power series solution in r _1 near infinity is very similar. 

For very special choices of the parameters one expects to find solutions of EYM-equations correesponding to 
gauge groups that are subgroups of G. In one simple case this is easy to see. 

Since every compact semi-simple (non-Abelian) Lie group has 5/7(2) as a subgroup there must be SU(2)- 
solutions embedded among general G-solutions for any symmetry group action. There may be several 
conjugacy classes of such subgroups, but there is a distinguished one related to the homomorphism defined 
by the action. Let the action and hence A be given and pick any f2 + thus selecting a specific st(2)-subalgebra 
generated by the triple {A , fi+, ft- = — c(il + )} among the conjugacy class associated with A . If we then 
let 

A+(r) = v(r)0+ with v(0) = 1 (3.24) 



it follows from (2.14) that v(r) = e no u(r) for a constant 70 and a real function u(r). The remaining equations 
fl2.llD - Q2.13l ) then become 

m' = gt(Nu' 2 + ±r- 2 (l-u 2 ) 2 ), (3-25) 

S^S' = 2g 2 r- 1 u' 2 , (3.26) 

r 2 Nu" + (2m - g^r" 1 ^ - u 2 ) 2 )u' + (1 - u 2 )u = 0, (3.27) 

where go := i||Ao||. They reduce with 

P ■= r/go, p := m/g (3.28) 

to the equations for the SU (2) theory, 

dpi/dp = N{du/dp) 2 + \p- 2 (\ - u 2 ) 2 , (3.29) 

S- l dS/dp = 2r-\du/dp) 2 , (3.30) 

p 2 Nd 2 u/dp 2 + (2fi - p- l {l - u 2 ) 2 )du/dp + (1 - u 2 )u = 0, (3.31) 

where now N = 1 — 2p,/p. 

Much is known about the solutions of these equations. In particular, it follows, that for any compact 
semi-simple Lie group G and any action of the symmetry group SU(2) an infinite discrete set of global 
asymptotically flat (soliton and black hole) solutions exists. 
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4 The spherically symmetric static EYM models for some small gauge groups 

We have not yet found a simple general way to derive the held equations for arbitrary semi-simple g and 
arbitrary choice of Ao. But in the following tables [| to || we list for lower dimensional Lie groups and the 
actions of SU (2) given by their characteristic some basic properties of the system of equations, namely the 
size of the size M = \S^ \ of the Stiefel set which corresponds to the number of complex functions w a (r)) 
that describe the gauge potential, the set £ with the superscript denoting the dimension of the eigenspace if 
it is greater than one, and the subalgebra to which the equations reduce if S\ is a II-systcm ('-' indicates 
that there is no reduction). We leave out the principal action, which is always regular, as well as the trivial 
action of the symmetry group. 

Just to give an idea of the structure the M-dimensional real subvariety S of C M we list in table ^ the 
equations that define it for a few cases. (In the regular case for a Lie algebra of rank I the equations would 
require that \w\\ = const for certain fixed constants for all i = 1, 2, . . . , t.) 

When one wishes to find global numerical solutions one can, for example, pick a simple choice for A + (oo) - 
it seems possible to choose all the w a {oo) real, and many of them 0. But then the data at r = must be 
left arbitrary, i.e. a suitable parametrization of the variety £ must be introduced. This is not difficult for 
the low dimensional cases, but not straightforward. 



Table 1: Irregular actions for Lie algebras A 2 to A 4 and B 2 to B4 



X 


l-?Aol 


£ 


reduction 


X 


\Sa \ 


£ 


reduction 




A 2 


(SU(3)) 


B 2 (50(5)) 


11 


1 


1 


At 


01 


1 


1 


A x 










20 


3 


l 3 






A 3 


(SU(4)) 


B 3 (SO(7)) 


101 


1 


1 


Ax 


010 


1 


1 


Ax 


020 


4 


l 4 




101 


2 


1=* 


Ax © Ax 


202 


4 


1 3 ,6 




220 


4 


1,2*, 3 












200 


5 


l b 












020 


6 


1 5 ,2 






A 4 


(SU(5)) 


B 4 (SO(9)) 


1001 


1 


1 


A 2 


0100 


1 


1 


Ax 


2112 


3 


1,2,3 


A 3 


2101 


3 


1 2 ,3 


B 2 ®Ax 


1111 


3 


l a ,2 


a 3 © a 2 


1010 


4 


l 4 




0110 


4 


2 4 




0201 


5 


1,2 3 ,3 




2002 


6 


2 b ,3 




2220 


5 


l,3 a ,5 












0001 


6 


1« 












2200 


6 


1,2 4 ,3 












2000 


7 


l 7 












2020 


8 


1 4 ,2 2 ,3 2 












0020 


9 














0200 


10 


l y ,2 
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Table 2: Irregular actions for the Lie algebras C 3 (Sp(6)), C 4 (Sp(8)) and D 4 {SO (8)) 



X 


ISaoI 


£ reduction 


X \Sa \ £ reduction 


C 3 (Sp(6)) 




100 


1 


1 












210 


2 


1,3 


c 2 










010 


3 


l y 












020 


4 


1,2 3 












202 


5 


1 3 ,2,3 












002 


6 


1« 












C 4 {Sp{8)) 


£>4 {SO{8)) 


1000 


1 


1 




0100 


1 


1 


Ai 


2100 


2 


1,3 


c 2 


1011 


3 


l 3 




0100 


3 


l y 




0202 


5 


1,2 3 ,3 




2210 


3 


1,3,5 


c 3 


0220 


5 


1,2 3 ,3 




2010 


5 


1 3 ,2,3 




2200 


5 


1,2 s , 3 




0110 


5 


1 2 ,2 3 




0002 


6 


l ti 




0010 


6 


l 6 




0020 


6 


1« 




2202 


6 


1 2 ,2,3 2 ,5 




2000 


6 


1 B 




0202 


7 


l 3 ,2,3 a 




2022 


6 


1 3 ,2,3 2 




0200 


8 


1 5 ,2 3 




0200 


8 


1 V ,2 




2002 


9 


l b ,2 2 ,3 












0002 


10 


!l0 













Table 3: Irregular actions for the Lie algebras F4 and G 2 



X 




C 

G 


reduction 


X 


\Sa \ 


£ 


reduction 




G 2 


1000 


1 


1 


A! 


01 


1 


1 


A! 


1012 


3 


1,3,5 


c 3 


10 


1 


1 


A x 


1010 


5 


1 3 ,2,3 




02 


4 


1 3 ,2 




0101 


5 


1 2 ,2 3 












0100 


6 


1« 












2001 


6 


1,2*. 3 












2202 


6 


1,2,3,5^,7 












0001 


7 


l 7 












2200 


7 


1,3 & ,5 












0002 


8 


1,2 V 












0202 


8 


1 3 ,2,3,4,5 2 












0010 


9 


1 6 ,2 3 












0200 


12 


l b ,2 4 ,3 2 












2000 


14 


l i3 ,2 
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Table 4: The "initial value" surface £ for a few irregular actions for some Lie algebras 



Lie algebra 


X 






090 


\wi\ 2 + \w 2 \ 2 = 1, \w 3 \ = \w 2 \, \Wi\ = \wi\, W1W3 = -W 2 U>4, 

W\W 2 = —W3W4 




1010 


\wi\' 2 + \w 2 \ 2 = 1, \w 3 \ = \wi\, \wi\ = 1, w 2 w 3 = wiw 2 


E e 


120001 


\wi\ 2 + \w 3 \' 2 + \w 5 \' 2 = 3, \w 4 \ = \w 3 \, \w 5 \ = \w 2 \, \w 6 \ = 

2, | W^j = \lVx\, U>iW 2 = W5W7, U>iW 3 = —W4W7, W1W4 = 
-W 3 W 7 , WiW 5 = W 2 W 7 , W 2 W 3 = w 2 w 4 = -w 3 w 5 


Fi 


1010 


\w!\ 2 + \w 3 \ 2 = 3, \w 3 \ 2 + \w 4 \ 2 + \w 5 \ 2 = 4, -| W2 p+2]™ 3 P+ 

1 W5 | 2 = 3, WlW 3 + W 2 W 4 — W4W5 


G 2 


02 


\ Wl \ z + 2\w 2 \ 2 + \wtf = 2, \w 2 \ A + 2\w 3 \ 2 + \w 4 \ 2 = 2, 
W1W2 + 2w 2 w 3 + w 3 w 4 = 



5 Algebraic Results 

In this section we collect all of the algebraic results needed to prove the local existence theorems. We will 
employ the same notation as in C3] section 6. 

Before proceeding, we recall some results from . 

Proposition 1. There exists M highest weight vectors S, 1 ,^, 2 , ■ ■ ■ ,£, M for the adjoint representation of q on 
g that satisfy 

(i) the £ J have weights 2kj where j = 1, 2, . . . , M and 1 = k\ < k 2 < ■ ■ ■ < kM> 

(ii) ifV{^) denotes the irreducible submodule of q generated by £ J , then the sum ^(£"0 * s direct, 

(iii) iftf = (l/l\)nl.& then 

<#?') = (-!)'&,-!. (5-1) 

(iv) M — \S.\\ and the set {^ 1 k ._ 1 \ j = 1,2, . . . M} forms a basis for V 2 over C. 
According to Lemma 1 of jl8| the R- linear operator A : g — > g defined by ( ft- 14 ) satisfies 

A(V 2 ) C V 2 (5.2) 

and hence restricts to an operator on V 2 which we denote by A 2 . 
We label the integers kj from proposition [j] as follows 

1 = fe/i = = • ■ ■ = fcjj+mj-i < kj 2 = kj. 2+ i = • ■ ■ = kj 2 + m2 -i 

<■■■ <k Jl = k Jl+ i = ■■■ = kj I+mi -i , 

where J\ = 1, J; + m; = J1+1 for I = 1,2, ... ,1 and J7+1 = M — 1. Then 

f = {k, :=k, h \l = \,2,... ,1}. 
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The set {£fe._i | j = 1,2, .. . M} forms a basis over C of V% by proposition [l] (iv) while I is the number of 
distinct nonzero eigenvalues of A 2 . Therefore the set of vectors {X l s , Y% |i = l,2,... s = 0, 1, . . . ,m;— 1} 
where 

X[ := { ?g£ ^ " ° dd and Y* := iX[ , (5.3) 

I Kk, -1 " 1S even 

forms a basis of V2 over R. Then lemma 2 of |l8| shows that this basis is an eigenbasis of V 2 and we have 
A 2 (X l s ) = k l {k l + l)X l s and A 2 {Y l s ) = for I = 1, 2, . . . , I s = 0, 1, . . . , m t - 1 . (5.4) 

An immediate consequence of this result is that spec(A 2 ) = {0} U {kj(kj + 1) | j = 1, 2, . . .7} and mj is the 
dimension of the eigenspace corresponding to the eigenvalue kj (kj + 1). Note that / is the number of distinct 
positive eigenvalues of A 2 . 
Define 

E l = span R {Y^ \ s = 0, 1, . . . , mi - 1} , E l + = sp&n R {X l s | s = 0, 1, . . . , mi - 1} , (5.5) 

and 

E = 04, E+=@E l + . (5.6) 

Then 

E =ker{A 2 ) (5.7) 



and 4 is the eigenspace of A 2 corresponding to the eigenvalue k;(k; + 1). Moreover, using proposition [l] 
(iv), it is clear that 

V 2 = E © E+ . (5.8) 
To simplify notation in what follows, we introduce one more quantity 

E l :=®l%®E* + . (5.9) 

q=0 

We then have the useful lemma from fl§| l. 

Lemma 2. If X E V 2 then X E E l if and only if ffy. X = or fl k ^ +2 .c(X) = 0. 
We will also need the the map ~: Z>_i — * {1, 2, . . . ,1} defined by 

~\ = 6 = 1 and s = max{ I | k ; < s } if s > 0. (5.10) 

It is shown in lemma 5 of [H that this map satisfies 

kg < s for every s E Z> and kjp < s < ks + i for every s E {0, 1, ... k/ — 1}. (5-11) 

The last result we will need from |p8j is the following lemma. 

Lemma 3. If X E V 2 , k p + s < k p+1 (s > 0), and Q,^ +S .X = 0, then Q, k £.X = 0. 
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We will also frequently use the following fact 

(l±l)~=l±l Ml EE. (5.12) 

Proposition 2. If X a G -E 5 , F 6 G B B , and Z c G £ £ i/ien [[c(X a ), F 6 ], Z c ] G £( a + b + c )~. 

Proof. Suppose X a E E d , Y b E E b , and Z c G E~ c . Then 

0^ +2 .c(X o ) = Vt^.Y a = fl^.Zc = (5.13) 

by lemma 0. Now, 



W + .[[c{X a ),Y b ],Z c ]=Y J Y.( P l ) (L) W > 

1=0 m=o ^ ' \ / 



where W pabcLm = \[Q.™.c(X a ), tf+ m .Y b ], S%T l .Z e ]. It then follows from ( p3| ) that ^ = if to > kg + 2 
or I — m > kg or p — I > kg. Thus W pa bclm = unless p < Z + kg < to + kg + kg < kg + kg + kg + 2. 
But this can never be satisfied if p = kg + kg + kg and so we arrive at fi + a b " \[c(X a ), Y b ], Z c ] = 0. But 
kg + kg + kg < a + 6 + c by ( |5.11[ ) and hence it follows that Q a ^ b+c .[[c(X a ) 1 Y b ], Z c ] = 0. But then lemma | 
implies that n^ a+b+cT .[[c(X a ), Y b ],Z c ] = and hence [[c(X a ), Y b ],Z c ] E £( a + b + c r. □ 

Proposition 3. If X a G E" andY b G £7 S ften [[X OJ c(Y b )], £}+], X a ],Y b ], [[fi_, JC ], Y 6 ] E E^' . 
Proof. This proposition can be proved using the same techniques as proposition |[ □ 

Lemma 4. If I E E and Z E E\® E 1 ' 1 then Q 1 ^ 1 .c(Z) = 1(1 + l)^ 1 .Z . 

Proof. Since Z G E l + © E l ~ x there exist real constants a*, such that 

l-l ro g — 1 mf— 1 

2=EEK+ <OT + E a f4 ( 5 - 14 ) 

q=l s=Q s=0 

where X* = ^-i if k 9 is odd and Xf = if k 9 is even. But 

Sl l + l .X* = for q< l-l (5.15) 

by lemma and so we get 

mj— 1 

fi^ 1 .Z= ^ afO^.xf. (5.16) 

s=0 

Now, c(XJ') = Cifj+i if kg is odd and c(Xj) = iC^+i k 9 is even by proposition |l| and so 

{l 2 + .c(XI) = k q (k q + l)X*. (5.17) 
Since I E E implies that kg = I, it follows easily form ( (5.14 ), (5.15), and ( |5.17 ) that 

m r -i 

n l +\ c {z) = i(i + a p l + 1 - x l ■ ( 5 - 18 ) 

s=0 

Comparing (|5.16[) and (|5.18| ), we see that fi'+ 1 .c(Z) = 1(1 + l)Q, l + 1 .Z . □ 
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Proposition 4. If Z\ with I — 0, 1, 



fc is a sequence of vectors such that 



Z a = n + , ZiE E> l = l,2,...k and Z t G E, 



ifle£ 



then for any j = 1, . . . , I 



k-i 



Err tr, \ rr 1 ^ 1 ^ if k 4 £ 

[W-),^],n + ]6| # _ 1 ^ 



Proof. Suppose is as in the hypotheses of the proposition, then 



kr+2 



c (Zj) = ni r .Zj = o . 



Now, 



fc-i 



fc-ip+i 

■fiL.3[^ i -,,z.] 1 n + ] = EE 

8=1 8=1 i = 



P+l 
Z 



Wpkls 



(5.19) 



(5.20) 



where WpW, = [O^ .c(Z fc -«), n£ +1 ~'.Z a ]. From ( pg| ) we see that W p ki s = if I > k(fc— «)- or p + 1 — / > kg. 
Thus 



Wpfc/ S = unless p+l<l + k s <k§ + k( fc _ s )~ + 2 



(5.21) 



Now, suppose p — fc. Then using the fact that kg < s and k( k _ s y < k — s, we get from fl5.2l| ) that W p fc; s = 
unless fc + l<Z + kg<fc + 2. Since this is impossible to satisfy W p ki s = for all /, s. Thus the sum (5.20) 
vanishes, i.e. — 12* . X)s=i [[ c (^fc-s> ^s]j = 0, and we get 



fe-l 



5^[[c(z fc _„zj,n + ]G^ A 



(5.22) 



8=1 



by lemmas || and |[ Suppose further that k £ £ and let p = fc — 1. Then a p fe; s = unless fc < I + kg < 
+ + 2 by ( |5.21 ). Now, kg < s and k(fc_ s )- < k — s, so suppose kg < s or k( fc _ s )~ < k — s. Then 

k(fc_ s )~ + kg < fc — s + s = k which will make the inequality k < I + kg < kg + k( k _ s y + 2 impossible to 
satisfy. Therefore we see that W p kis — unless k( fc _ s )~ = k— s and kg = s (i.e. k — s,s G £). However, if 
k(fc- s )~ = k — s and kg = s, then k < I + kg < kg + k(/._ s v- + 2 will satisfied only if Z + s = 1. So W p hU = 
unless k — s, s € £ and Z + s = 1. This allows us to write the sum (5.20) as 



fc— l 



-n*- 1 .5^[[ c (z fc _ i ),zj,n + ] 



fc-i 

E 

8=1 
fc-1 



E 



fc - s + 1 
fc! 



- (k- 8-1)1(8-1) 



(fc - s)(fc - a + l)^-'- 1 .^.., Q'+ l .Z a ] 



by lemma 
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Assume now that k is odd. Then we can write the above sum as 

fc-i fc-i 



-0^.]r[[c(Z fe _ s ),Z s ],0 + ]= £ ___«__[ n ^i.z^ - ,^.Z.] 

fe— 1 

2 fc 1 

+ E (A _ s - i)! (s _ i)! «+ 
= E ( fc - a -i)!( s -i)! {[^ 1 -^-^r 1 -^] + [nr 1 -^.n+"*" 1 -^-]} = • 

Similar arguments show that — fi* -1 . J2^Zi[[ c {Zk-s, Z s ], fi+] = if fc is even. Therefore Es=i[[ c (^i-s)' -^h ^+] e 
J B fc - 1 by lemmas | and | and G5.12| ). □ 

Proposition 5. J/ wii/i I — 0,1, ... k is a sequence of vectors such that 

Z a = n + , Zi 6 E l l = l,2,...k and Zi E E l + E 1 ^ 1 if I € £ 
then, for any j = 1, . . . ,1 

V^rr (7 ^ 71 7 1 c J # ^ £ 

E E^-^ ^Ufe-j] e i pit-i ... 



Proof. Proved using similar arguments as for proposition |} □ 

6 Local Existence Proofs 

For q = 1,2,... I let 

pr q + : V 2 ^ E q + , pr q :V 2 ^E q , and pr ? : V 2 - E q © B« (6.1) 
denote the projections determined by the decomposition (pTq), (t5.6|) of V%. 



6.1 Solutions regular at the origin 

Theorem 2. Fia: X E E + and fl + G £/ia£ satisfies [fL|_,0_] = Ao where f2_ := — c(fi+). Tften i/iere 
ea;«st a unique solution { A_|_ (r, Y), m(r, Y")} to £/ie system of differential equations (2.11) and ( 2.1 3| ) £/wrf is 
analytic in a neighborhood of (r, Y) = (0, X) iralx £7+ cmd satisfies m = 0(r 3 ) and 

p4(A+ - = F s r s+1 + 0(r s+2 ) , prg(A+) = 0{r s+2 ) Vs e £ 

where f2+ := A+(0) and Y s :— pr5_(Y). Moreover, these solutions also satisfy P = 0(r 4 ) and G — 0(r 2 ). 

Proof. Introduce new variables { u+ , u® | s & £ } via 

uf :=pr»(A+-n+)r— 1 and u° s := pr^A+K*- 2 (6.2) 

where J7 + := A+(0). This allows us to write A + as 

A+(r) =n + + 5>+(r) + ru2(r)) r s+1 . (6.3) 

But regulari ty at r = r equire s that f2+ satisfy [f2+,f2_] = Ao where f2_ := — c(fl+). This can be seen 
easily from ( 2.16| ), ( 2.24 ), ( 2.25| ) and the requirement that the pressure remains finite at r = 0. 
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Lemma 5. For every s <= £ the re exists analytic maps : E + — > -Eg and F 2 : E x E + xl^ Eq(BE + 



such that, for T given in ( 2.17 ), \>y s T = —s(s + l)ufr s+1 + r s+2 J 7 }(u + ) + r s+3 T 2 (u° , u + , r) where 



i° -=^2 u a and u+ : =^2 u t 



(6.4) 



Proof. Let u s = rvP s + u+. Then from ( 2.17 ) we find 



T = Y J Mu a )r a+1 + \Y, ([[u a ,c(u b )},n + ] + [[a + ,c(u a )},u b ] + lin-,u a ],u b })r a 



+b+2 



a££ 



a,b££ 
1 

+ 2 



\ [[ua,c{u b )lu c ]r a+b+C+ * 



But from (p]j) we see that 



M{u a ) = a(a + l)u+ 



+ 



Also, for a & £ we have ka = a by lemma 5.11, and hence 

a < b a < b . 



(6.5) 



(6.6) 



Using the two results ( p.5|) and (6J3), we get 



pr 



~<T = -s(s + l)u+r s+1 + IY1 pl '" ( c ( u »)]. 0+] + [[«+, c(« a )],«6] 



-[[ri_,M Q ],u b ] r 



„a+b+2 



a+fe+c>s 



[u a ,c(u b )],u c ] )r 



a+b+c+3 



by propositions || and |^. Substituting u a = rvP a + u+ into the above expression completes the proof. □ 
For every s G £ define 

v+:=uf and v° s := (ru° s )' . (6.7) 



Lemma 6. There exists analytic functions P : Eq x E + xl^l and G : E x E x _E + x £7 + x 
sucft i^P = r 4 || M +|| + r 5 P(u°, u + , r) and G = r 2 u^" || + r 3 G(u°, i>°, it + , t> + , r) where 

V° := ^ v a and V+ : = V t ■ 



(6.8) 



and u° , u + are defined by (6.4). 



Proo f. Th e ex istence of the analytic function G follows easily from the definition ( 2.15 ) of G and equations 
~~ | and (HJ). 



Similarly from the definition (2.15) of P we have P = ^-JjT2 + , c(u^)] + [il_,uf]\\ +r 5 Q(u°,u + ,r) where Q 
is a polynomial in r, u° and u + . Now, using (2.23), (3.14), and ^(it^) = 2u± , it is not difficult to show 



that | I [n + , c(ttj")] + = and this completes the proof. □ 
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From (5.5)-(5.7) we get A^jut) = s(s + and A 2 (u° s ) — sinc e u+ 6 El and u° € Eft. Using this result, 
lemma |5|, and equations ( |6.4| ) and ([T^), the field equations ( [2.1l| ) and ( 2.13 ) can be written as 



ruj = rv„ 



= -2(s + l)v+ - — ( m --p)v+- 



s(s+ 1) / 1 



2(5+1) 

r 2 N 



N- 1 



_-l)pr»^(«+)-pr»^(« + ), 



: -2(s + 1K° - s(s + l)u° - 1 (m - ~p) ^ 
2(s + 1) 



riV 



m--P)«2--pi^*(u ,u+ r) 



--l) pr^ 1 ( W +)-pr^ 1 ( W +), 

where s S £. For every s S £, introduce two new variables 

z s := -( s + l) u o_i±l v o 

and define 

1 



(s + l)u° s + v° s , 



fs 



2 

"iV 



m P I 

r 



riV V r 



1 



-prg^VV) - ( - - 1 ) prgj*(u+) 



AT 



Then equations (6.11) and (6.12) can be written as 

rx' s = -(s + 2)x s + £±lprg^(t i +) + 



s + 1 



-( S + 2) 2/s -pr^ s 1 (u+)/ S) 



for every s € E . Define fi = 



m — r 3 ||ux"|| j- Then the mass equation ( 2.11 ) becomes 

rp! = -3/i + r |-P(u°, u + , r) + G(u°, v°, u + ,v + , r) - 2(( u+\v+)) 

-2r (2|| M +|| 2 +rG(M ,i> ,u + ,w + ,r)) } . 



For every s € £, introduce one last change of variables 



v' := i> 



2(s + 1) 



pr s Q J r l(u + ) , x s 



1 



s(s + 2) 



proJ r s 1 (u + ), and y s :=y s 



(6.9) 



(6.10) 
(6.11) 



(6.12) 



(6.13) 
(6.14) 



(6.15) 



+ 



Let w 



+ — 



see "« ' 



= £ se£ £ s , y : = Eaeef* and ^( r ) : = (x(r),y(r),u+(r),v + (r),n(r),r). Fix X e E + 
and let Mx be a neighborhood of X in P + . Define a set D(Nx, e) by £ )( A/"x, e ) : — Pp x P x ATy x P+ x 
(— e, e) x (— e, e). Then using lemmas || and || and equations ( |6.9[ ), ( |6.1C| ), ( 6.1 3| ) , ( |6.14 ), and (|6 .15 ) one can 
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show that there exists an e > and analytic maps U S ,V S ■ D(Afx, e) —* El , X s ,y a : D(Afx, e) — ^ -Eo > an d 
M : D(Af x ,e) -> E such that for all s e £ 

ruf\r) — rU s (r](r)), rvf\r) — — 2(s + l)t)^(r) + rV s (7y(r)), 
rx a '(r) = -(a + 2)x s (r) + rX a (r)(r)), ry s '{r) = -(s + l)i s (r) + r^^r)), 

and r/i'(r) = — 3/x(r) + rA1(??(r)). This system of differential equations is in the form to which theorem [l] 
applies. Applying this theorem shows that for fixed X € E + there exist a unique solution {uf(r, Y),vf(r, Y), 
x s (r,Y),v s (r,Y), fi(r,Y)} that is analytic in a neighborhood of (r,Y) — (0,X) and that satisfies uf(r,Y) = 
Y s + 0(r), v+(r,Y) = O(r), x s (r,Y) = O(r), y s (r,F) = O(r), and /x(s,F) = O(r) where Y" s = p4(F). From 
these results it is not hard to verify that m(r) = 0(r 3 ) and u°(r) = 0(r°) Also, it is clear that P = 0(r 4 ) 
and G = 0(r 2 ) by lemma |[ □ 



Theorem 3. Every solution from theorem^ satisfies equation ( ^.14 ) on a neighborhood of r = 0. 



Proof. Let {A + (r), m(r)} be a solution of the equations ( 2.11 ) and (|2.13|) on a neighborhood A/" of r = 0, 



which we know exists by theorem H. From theorem [2|it is clear that A' + (0) = 0, and this implies that 7(0) = 
where 7(7-) is defined in lemma [l|7 Also, because m = 0(r 3 ) and P — 0(r 4 ) for these solutions we see, by 
shrinking Af if necessary, that the function f(r) = —2(r 2 N)~ 1 (m — ^P) is analytic on Af. From lemma 
7 satisfies the differential equation 7' = f{r)^f. Solving this equation we find j(r) = 7(0) exp(J^ f(T)dr) for 
all r £ Af . But 7(0) = 0, so j(r) = for all r € Af. □ 



6.2 Asymptotically flat solutions 



In proving that local solutions exis t nea r r = 0, w e were able to "guess" the appropriate transformations 
needed to bring the field equations ( 2.11 ) and ( 2.13 ) in to a form for which theorem [l] applies. Near, r — 00 
the equations become much more difficult to analyze and guessi ng th e appr opria te transformation is no 
longer possible. Instead, we will show that the field equations ( 2.11 ) and ( 2.1 3|) admit a formal power 
series solutions about the point r = 00 . Th is for mal p ower series solution will then be used to construct a 
transformation to bring the equations ( 2.11 ) and ( 2.1 3| ) in to a form for which theorem [l| applies. 



Let 



and define 



/ = - 

dz 



for any function /. Then the field equations ( 2.11 ) and ( 2.13 ) can be written as 

NG + z 2 P = 0, 



1 
z m 



z 2 NA + + 2z(l - 3zm + z 2 P)A + + T = . 



Assume a powerseries expansion of the form 

00 

A+ = J2 A +*z 

We will define A_.fc := — c(A 
vanish we have that [f2+,f2_] = 
yields the recurrence equations 



and 



fc=0 



^2m k z k 

k=0 



(6.16) 
(6.17) 

(6.18) 



if) and fl± :— A±^. From the requirement that the total magnetic charge 
= Aq. Substituting the powerseries ( 3.18 ) in the equations ( 6.16 ) and ( 6.17 ) 



mi = TO2 = , TOfc = — 

k 



1 I \ 

Y -Gfc-3 + 2) mjGk-j-4 - Pk-i 

V / 



k = 3,4,5, 



A 2 (A + , fc ) - k(k + l)A+, h = h k + f k k = 1, 2, 3, 



(6.19) 
(6.20) 
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where 



k 

G k ■= \ X>' + !)( fc + 1 - A+, fc+ i-i|A +)j+1 )) k > , 

k j 

F :=0, F fe :=-^^[[A_ J _ s ,A + , s ],A + , fc _ J ] fc > 1 , 

j=l s=0 

fc-1 

P = Pi=0, P fc :=-^((^|F fe _,)> fc>2, 



A— 1 



^i=0, /i fc :=2^(fc-i-l)(P i _2-(fc-i + l)m i )A +>fe _ j _i k>2. 



(6.21) 
(6.22) 
(6.23) 
(6.24) 



and 



fc-1 j 



k-i 

v . fe - J ']+^[[A_, fe _ s ,A + s],n+] [> fc>2 



(6.25) 



^ j=l s=0 

Note that with these definitions that F — J2kLo ^kZ k and P = J2kLo PkZ k while G = J2kLo G , / £ z fe+4 . 



Theorem 4. Fix X £ E + andm^ £ R. TTien t/iere exists a unique solution {A +j fc, TOfe}^L to the recurrence 
equations (6.19) and (6.20) that satisfies 



m = moo , mi — m 2 — ; pr + A +> fc = X k V k £ £ 



and 



A +,fc e 



# i/ fc ^ £ 



(6.26) 



(6.27) 



where X k ■= pr5_A". 



Proof. Fix AT £ £7+, £ M, and let AT^ = pr*^A" for all k £ £. We will use induction to prove that the 
recurrence equations ( |6.19 ) and (3.2C) can be solved. When k = 1, the equations (6.19) and (6.20) reduce 
to mi = and A2(A +! i) — 2A +! i = 0. This can be solved in Eq © E\ by letting A+ ) i = X\. Note that since 
minf = 1, we have 1 = 1. 

We now assume that for k < I , {A+^mfc} is a solution to the recurrence equations ( |6 . 19 ) and (6.20) that 
satisfies 



A+,fc e 



E k 



if ^ £ 



E\ © E k ~ x if k £ £ 



It is clear from ( |6.20 ) that m; + i is then determined. From (6.21)- fl6.25|) and propositions || and || it follows 
that 



+ fi+i € 



P [+1 if Z + 1 ^ £ 
E f if Z + 1 e £ 



(6.28) 
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Equation ( 6.20 ) implies that 

(A 2 - (I + 1)(J + 2)1 y 2 )A+,, +1 = h l+1 + f l+1 . 
Suppose l + l<££. Then A 2 - (I + + 2)1 is invertible and 

= (a 2 + m + 2)1 y\h l+1 + f i+1 ) . 



(6.29) 



But then ( |6.28[ ) implies that A + ,;+i G 

Alternatively, suppose 1 + leS. Then ker(A 2 - (I + I) (I + 2)1 ) = E^ 1 by (fj)-|[7|) and dUl)- Therefore, 
d6.28| ) shows that 

A+, i+ i = ((A 2 - (/ + + 2)1 )| £ f) _1 + / i+x ) + 

solves (|6.29D since AT i+ i G ker(A 2 - (I + + 2)1 ). It also clear from and S that 

A+,(+i G -Eo ® -^+- This prove that A+^+i satisfies the induction hypothesis and so the proof is 

complete. □ 

Theorem 5. Fix X G E + , > 0, and f2 + 6 i/iaf satisfies [fi + ,fl_] = Ao where f2_- := — c(f2+). 
TTien i/iere eiisf a unique solution {A + (r, a, Y~), m(r, a, K)} to £/ie system of differential equations (2.11) and 
( 2.1 3j ) £/ia£ is analytic in the variables (r _1 ,a, Y) in a neighborhood of (0,rrioo, X) G K 2 x E+ and satisfies 
m = a + 0(r~ 3 ) and 

Y, 



pr;(A + -O + ) = ^ + 0(r- s ), prg(A + ) - 0(r- s ) Vs G 5 



where 



+ 



A+(0) and Y, := pri(Y~). 



Proof. Fix A" G E+ and let A +J t = A + _k(X,m 00 ) and = mk{X, TOqo) be solutions to the recurrence 
equations (6.19) and (6.20) which satisfy mo = TOoo , TOi = m 2 = 0, pr^A +J t = for all k £ £, and 

A+.A; G 



if k i £ 



k 

q=l 

E\ ® ® k ~\ El © El iffcGf 



where Xk = pr+ AT. Define U :— J2k=o A+,fc^ fc and M := Y^k=o m k zk an d introduce new variables 4>(z) and 
<t{z) via 

A+ = U + z n - 3 cj) 
where the integer n is to be chosen later. Define 



and m = M + z n 3 a 



1 - 2Mz , F p 



UA + [U,c(U)]), 



F p :=-i[F p ,U], P p := \\\Fh p \ 2 , and G p := ±z 4 \\U\ 



(6.30) 

(6.31) 
(6.32) 



From these definitions it is clear that the quantities N p , F p , T Pl P p and G p are all polynomial in the variables 
X and moo- Now, because U and M are the first n terms in the powerseries solution to the field equations 
fl6.16| ) and ( 6.17| ) about the point z = they satisfy 

z 2 N°U + 2z(l - 3zM + z 2 P p )U A T p = z n - 1 {a 1 {X, m^) + a 2 [X,m 00 )y) , 
z 2 M + N P G P + z 2 P p = z n b(X, moo) , 
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where a±, a% : V 2 x 



and b : V 2 x 



are polynomial in their variables. 



From ( |6.32| ) and ( |2.17D it follows that T = T v - z n ~ 3 A 2 ((j)) 



Y^j=i^R,j(4>,X,m !X ,z) where T R , 



V 2 x E + xRxR^R j = 1, 2, 3 are analytic maps that satisfy J-Rj(eYi,Y 2 , %i, x 2 ) = e° Tr ,j(Yi, Y 2 , xi,x 2 ) 
for all e € IR. It is also not difficult to see from ( |6.32D and ( |2.17D that 



G = Z — ||C/|| 2 + 2 4 «L>|(n - 3)z"-> + z"~ 3 0)) + ~||(n - 3)z"- 4 + z n ~\ 



and P = P p + z n 2 PR,j((j>, X, m^, z) where P R<j : V 2 x E^ 



1, 2, 3, 4 are analytic 



maps that satisfy Prj(cYi,Y 2 , xi, x 2 ) — Prj(Yi 7 Y 2i x 1: x 2 ) for all eel. Note also that N = N p — 2z n 2 a. 
Let 



uj := (f) and 9 := z \ 



(6.33) 



Using the above results, straightforward calculation shows that there exists anal ytic maps Q : V 2 x V 2 x E+ x 
M. 3 — > V2 and M : V 2 x V 2 x E + x M 3 — > R such that equation (6.16) and (6.17) can be written as 



and 



z<7 = — (ri — 3)cr + z 2 .M w, X, rrioo, cr, z) , 
Z# = -9 + lu , 

7V( Z £ + 2(n - 3)uj + (n- 3)(n - 4)0) + 2(n - 3)6» 

+ 2w - A 2 (0) - zg(9, w, X, moo, a, 2) = . 



We can rewrite (6.3E) as 

= -2(n- 2)w + (A 2 ~ (n- 3)(n- 2)1 )0 + z£(6>, w, X, m^, cr, 

where 



(6.34) 
(6.35) 



(6.36) 
(6.37) 



G(Yi, Y 2 ,Y 3 , xi, x 2 , x 3 ) 



x 3 \1 - 2(M(Y 3lXl ) + x^- 3 x 2 )x 3 

2Y 2 -A 2 (Y 1 ))+g((Y 1 ,Y 2 ,Y 3 ,x 1 ,x 2 ,x 3 ) 



1 I (2(n-3)n+ 



Because is analytic, it is clear that Q is analytic in a neighborhood of (0, 0, X, moo, 0, 0) £ V 2 xV 2 x E + x M 3 . 
For Y € V2 and s S £, define V s + := pr^_y and Y s ° := prgY. Recalling that prgA2 = A2PI0 = and 
pr^^4 2 = A 2 pr^ = s(s + l)prj_ for every s € £, we can write ( |6.35| ) and ( |6.37 ) as 



zw+ = ~2(n - 2)cj+ + (s(s + 1) - (n - 3)(n - 2))6»+ + zg+ s (6>, w, X, m^, cr, z) , 
z0? = -0?+ W ° 



and 



z£° = -2(n - 2)cj+ + -(n - 3)(n - 2)6>° + z£°(0, w, X, m^, cr, z) 
for all s £ £. For every s € £, introduce one last change of variables 

(I := (-n + 3)9° u ° , £ == (n - 2)0° + ^ , 



(6.38) 
(6.39) 
(6.40) 

(6.41) 



Vs ■ = 



1 

2s + 1 



((s-n + 3)0+ r? 2 



1 

2s + 1 



((a+n-2)0+ + w+) 
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and let ^ '-=^2 s e£ anc ^ ^ := V J S 3 — 1j 2. Using this transformation we see from the above results 

that ( |6. 38 - 6. 41 ) can be written as 



*Cl = - (« " 3 Ki + zlCl (C 1 , C 2 , r? 1 ,V 2 , X, m x , a, z) 

zr[> s = -(n- (-l) J s-j)rf 8 + zW s {C} , ( 2 , T) 1 , n 2 , X, m^, a, z) 



Vse£,j = l,2, (6.42) 
Vs6f ,j = l,2. (6.43) 



where KP S and 7i J s (j = 1,2) are i?o and valued maps, respectively, that are analytic in a neighborhood 
of (0, 0, X, moo, 0,0) E V 2 x V 2 x E + x M. 3 . The s ys tem o f differential equations given by Q6.34| ), ( |6.42| ), and 
( ]6.43D is equivalent to the original system ( |6.16j ), ( |6 . 17[ ) . Moreover, if we choose n — max{3, 3 + max£}, 
then fl6~34[ ), §J% ), and ( p5| ) are in a form to which theorem |l| applies. Applying this theorem shows 
that there exist a unique solution {o~(z,a,Y),Q(z,a,Y),rf s (z,a,Y)} that is analytic a neighborhood of 
(z,a,Y) = (0,moo,X) and that satisfies Q(z) — O(z) , ( 3 s (z) = O(z), and a(z) = O(z). It then follows from 



flOOD that p4(A+ - n+) = r s z s + 0(r s ), prg(A+) = 0(z s ), and m 



0(z 6 



□ 



Theorem 6. Every solution from theorem^ satisfies equation ( 2.14 ) on a neighborhood of r 1 = 0. 

Proof. Let z = l/r and {A + (z), m(z)} be a solution of the equations ( 2.11 ) and ( [2.13 ) on a neighborhood AT 
of z = 0, which we know exists by theorem |^. From lemma [l] it is easy to see that in terms of the z variable 



7 (z) = -z 2 ([A + (z), A_(z)] + [A_(*),A + (*)]) 



(6.44) 



and 7(2) satisfies 7(2;) = f(z)j(z) where f(z) — 2N 1 (m—2zP). By theorem || and shrinking Af if necessary, 
we see that f(z) is ana lytic on Af . Therefore we can solve ■y(z) = f(z)j(z) to get 7(2:) = 7(0) exp( /(r)dr) 
for all z G W . But ( |6.44 ) shows that 7(0) = 0. Therefore 7(2) = for all z G Af and the theorem is 
proved. □ 



6.3 Regular black hole solutions 

Theorem 7. Let t = r — ru and suppose X £ V 2 satisfies 



v{X) :=- -||A + [X, C (X)]|| 2 >0. 



Then there exist a unique solution {A + (i, Y), N(t, Y)} to the system of differential equations (2.11) and 
( 2.1 3j ) that is analytic in a neighborhood of (0, X) in K x V% and satisfies 

N(t) = v(Y)t + 0{t 2 ) and A+(t) =Y+ 0(t) . 

Proof. The proof is exactly the same as the proof of theorem 6 of with E + replaced by V?. □ 

Theorem 8. Every solution from theorem^ satisfies equation ( [2.13 ) on a neighborhood of r = ru . 

Proof. Let t = r — rn and suppose X € V% satisfies v := rn^ 1 — r^ 1 ||Ao + [X, c(X)]|| 2 > . Then we know 
by the previous theorem that there exists a solution {A+(t), N(t)} to the system of differential equations 
(2.11) and ( 2.13Q that is analytic in a neighborhood of t = and satisfies 

N(t) = vt + 0(t 2 ) and A+(t) = X + 0(t) . (6.45) 
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Figure 1: Gauge group 50(5) or B 2 : characteristic (20), A = 2h x + h 2 , (w a (0)) = (0,1,0), ||fi+|| 2 = 2, 
P = (0.1, —0.35, —0.28). The function v 2 is identically 0, 113 — —u\, and U3 = —vi. The quantity L is ||A + || 2 . 
For a globally regular solution it should nowhere exceed its value at the center and at infinity. 



Let f(t) := -2t((t + r H ) 2 N(t)) (m(t) - (t + r H y 1 P(t)) . Then (|6.45| ) shows that f(t) is analytic in a 
neighborhood of t = 0. A short calculation shows that /(0) = — 1, and therefore we can write f(t) = —l+tg(t) 
where g(t) is analytic in near t = 0. Consider the differential equation on V2, 

t^ = -r,(t)+tg(t)r 1 (t). (6.46) 

It has ?y(i) = as a solution, and therefore this is the unique analytic solution near t = by theorem |l|. But 
lemma |l] shows that -f(t) — [A+(i), ^jp(i)] + [A_(t), ^ t (t)] also solves ( 6.46| ) in an neighborhood of t = 0. 
Because 7(t) is analytic near t = 0, we must have j(t) = near t — 0. □ 



7 A numerical example 

So far we have not yet found a global numerical solution that has the correct fall off behavior for r — ► 00. 
The following example is for the gauge group SO(5) (B2) for the action with characteristic (20), one of the 
simplest irregular cases that will not reduce. Figure 1 shows a solution near the center with real initial 
values w Q (0) which develops nonzero imaginary components v a . Figure 2 shows a solution for large r. As 
is apparent the function decreases from its value 1 at infinity as r decreases down to some minimum, 
but then increases rapidly. Only by very careful tuning of the data at infinity one can possibly avoid this 
behavior and construct a 'physical' solution in which A^ alwas remains between and 1. For a globally 
bounded solution we have also indications that necessarily ||A + || < ||0+||. 
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-2-10 1 2 3 



4 



— v 2 

— v 3 



Figure 2: Gauge group 50(5) or B 2 : characteristic (20), A = 2hi + h 2 , («J«(oo)) = (0, 1,0), ||0+|| 2 = 2, 
Tiioo = 0.52, a = (—8, —2, —1). The function i>2 is identically 0, W3 = — u\, and W3 = — v\. The quantity L is 
||A + || 2 . For a globally regular solution it should nowhere exceed its value at the center and at infinity. Also 
N should remain between and 1. 
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